Semiclassical Theory of Chaotic Quantum Transport 
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We present a refined semiclassical approach to the Landauer conductance and Kubo conductivity of 
clean chaotic mesoscopic systems. We demonstrate for systems with uniformly hyperbolic dynamics 
that including off-diagonal contributions to double sums over classical paths gives a weak-localization 
correction in quantitative agreement with results from random matrix theory. We further discuss 
the magnetic field dependence. This semiclassical treatment accounts for current conservation. 
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Among the prominent wave phenomena which consti- 
tute mesoscopic behavior of small phase-coherent con- 
ductors, weak localization (WL) represents one key mech- 
anism. This quantum effect shows up as a decrease in 
the average conductivity with respect to the classical 
one. WL, originally observed for light and electron waves 
in disordered samples, has been extensively studied dur- 
ing the last decade for electrons in ballistic conductors, 
i.e. mesoscopic cavities or quantum dots where the elas- 
tic mean free path is considerably larger than the sys- 
tem size. WL is attributed to constructive interference 
of electron waves which are either coherently backscat- 
tered from impurities in disordered systems or multiply 
reflected at the boundaries of ballistic devices. 

In describing ballistic transport semiclassical methods 
have attracted much interest, since they establish a direct 
link between quantum transport properties and features 
of the corresponding classical dynamics, e.g. chaotic, in- 
tegrable, or mixed behaviour |^J^. This was demon- 
strated for clean microstructures in a seminal semiclas- 
sical approach Q to the average reflection in the Lan- 
dauer framework: the WL peak profile (as a function of 
a magnetic field) was shown to be Lorentzian for chaotic 
systems, probing in an impressive way the imprint of 
the classical dynamics on a measured quantum effect [D . 
However, while the line shape of the WL peak agreed 
with results from random matrix theory (RMT), the ap- 
proach turned out to be inadequate to give the correct 
WL magnitude for ballistic systems, contrary to the diffu- 
sive case 1^ . This long-lasting problem to semiclassically 
obtain the correct leading-order quantum correction to 
the conductance is related to the so-called diagonal ap- 
proximation used. It is based on the consideration of a 
restricted class of pairs of paths built from each backscat- 
tered orbit and its time-reversed partner, which moreover 
violates current conservation Similarly, WL is not 
captured in a related semiclassical approach to the Kubo 
conductivity of ballistic systems . 

The possible relevance of pairs of non-identical 
backscattered paths, differing slightly in their initial di- 
rections, was first pointed out by Argaman Later, 
Aleiner and Larkin j|] approached the problem of ballistic 



WL using both, perturbation theory and supersymmetri- 
cal methods to derive a RMT result for the conductance. 
However, their techniques still rely on the presence of 
quantum scatterers (to regularize the Liouville operator) 
and strictly speaking do not treat the case of a clean, 
disorder free, system. Their approach was semiclassically 
interpreted in Ref. arguing that diffraction-induced 
small angle scattering is relevant for ballistic WL. 

Here we present an adequate, current-conserving semi- 
classical treatment of the problem to quantitatively de- 
scribe the average quantum conductance in clean chaotic 
systems without relying on any diffraction or impurity 
scattering effects. We consider the leading-order off- 
diagonal contribution in a semiclassical loop expansion 
of the Landauer conductance. The relevant off-diagonal 
terms consist of pairs of orbits which are very close al- 
most everywhere (in configuration space), and differ only 
in whether they undergo or avoid a self-intersection with 
small crossing angle . Our results are strictly derived 
for chaotic systems with uniformly hyperbolic dynamics, 
but related results for ballistic cavities indicate |Q that 
they apply to general chaotic systems. 

We first compute semiclassical conductance contribu- 
tions beyond the diagonal approximation in the Landauer 
framework and later return to the corresponding prob- 
lem in the Kubo formalism. Consider a two-dimensional, 
classically chaotic clean cavity with two leads of width 
w[w') attached that support N{N') open channels. The 
Landauer formula for the conductance G then reads 
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Here tnm{E,B) denotes the transition amplitude be- 
tween incoming and outgoing channels m and n at en- 
ergy E in the presence of a magnetic field B. We first 
consider the case of time-reversal symmetry, B = 0, and 
return to the B-dependence of WL later. We assume 
that the ergodic time is much smaller than the escape 
time T of the cavity and that contributions from direct, 
lead-connecting processes are negligible. Then the fol- 
lowing RMT results for the transmission and reflection 
amplitudes hold which we give for later reference ■ 
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Our conductance calculation is based on the semiclas- 
sical representation of transmission amplitudes ||l8|] , 
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The sum runs over all lead-connecting trajectories 7 
which enter into the cavity at (a;,?;) with an angle 
sin 0m — fhiT I {kw) and exit the cavity at (x',y') with 
angle sin0,-i = fm / [kw'), where n = ±n, and p = hk is 
the momentum, see Fig. |l|(a). In Eq. (||), Sj is the clas- 
sical action, an element of the stability matrix, and 
= sgn(n)sgn(m) exp[z7r(m?//u> — riy'/w' — /i-y/2 + 1/4)] 
is a phase factor where contains the Morse index. An 
expression corresponding to Eq. (^) holds true for r„m in 
terms of paths reflected back. 

The Landauer Eq. contains products inmCm which 
semiclassically amounts to evaluate double sums over 
an infinite number of trajectory pairs. In a treatment 
of the energy-averaged conductance most pairs, con- 
sisting of orbits with uncorrelated actions, will cancel 
each other upon summation. The existing semiclassi- 
cal approach |^ is based on the diagonal approxima- 
tion, where only pairs of identical orbits or orbits related 
to each other by time inversion are taken into account. 
Then the phase factors from Eq. (jj) cancel, and one has 

l^nmlLg = W(2W)ET,(f.,™) |COS0'^COS0^MT'J-1. ^ 

First we give an alternative further evaluation of this 
expression employing the sum rule ^9 20 
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The sum runs over all orbits with periods T-y, which begin 
and end in intervals dy' and dy around y' and y with fixed 
orientations of the initial and final velocities. is the 

energy surface in phase space; S(-E) = 21:^1 A for billiards 
of area A. The factor p{T) - exp(-T/r) (for T ^ 00) 
accounts for the exponential loss of particles with velocity 
V which escape through the openings characterized by the 
escape rate 
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Upon applying the sum rule to the diagonal contri- 
bution, integrating over the lead cross sections, and in- 
cluding a factor 4 for each tupel (n, m) one finds for the 
transmission coefficient of an ergodic system 
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Correspondingly, the quantum reflection coefficient 
reads in the diagonal approximation 
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The semiclassical evaluation at this level yields the {k = 
0)-term of the RMT resuh (Hjl). Note that the use of the 
sum rule (|^) allows us to compute individual transmission 
and reflection coefficients, while Ref. [|j gives results only 
for the entire classical transmission and reflection. 

Summing the first term in Eq. (^) over all channels 
yields the classical reflection R = N'^/{N + N'). The se- 
cond term in Eq. (||) arises from contributions to |r„np 
from backscattered orbits paired with their time-reversed 
partners (elastic enhancement). This gives rise to the 
diagonal contribution to WL, 5i?diag = N/{N + N'), as 
derived in . In the limit N = N' 00 one has (5i?diag — 
1/2, deviating from the RMT resuh (5i?RMT = 1/4. 





FIG. 1. Sketch of an off-diagonal trajectory pair formed by 
a self-intersecting classical orbit (solid line) with small cross- 
ing angle e and a neighboring orbit (dashed) differing mainly 
in the region around the self-intersection. The paths represent 
orbits with many reffections at the system boundaries. They 
contribute to the quantum transmission (a) and reflection (b) . 

In the following we go beyond the diagonal approx- 
imation and consider pairs of different trajectories as 
sketched in Fig. Q in coordinate space. They consist of a 
long, self-intersecting orbit (solid Hne in Fig. |l|a),b)) with 
small crossing angle e forming a closed loop and a second 
orbit in its close vicinity (dashed line). The two orbits 
follow the two open trajectory segments, beginning and 
ending at (exponentially) close points at the lead mouths, 
in the same but the loop in opposite direction. Given a 
self-intersecting orbit with small e we showed that the 
neighbouring orbit indeed exists by linearizing the mo- 
tion in its vicinity |l^ . The action difference A5(e) of the 
two orbits is computed by expanding the action around 
the self-intersecting orbit up to second order. The re- 
sulting formula for AS', expressed through the elements 
of the stability matrices for the loop and the two open 
segments, is rather involved. Hence, we will focus from 
now on onto chaotic systems with uniformly hyperbolic 
dynamics characterized by a single Lyapunov exponent 
A and no conjugate points. One then finds [p^ 
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Since a partner orbit is associated with each self- 
intersection with smah crossing angle, we compute the 
conductance contribution from all such orbit pairs by 
first summing for each orbit "f(n,rh) over all e self- 
intersections and finally by integrating over e. Using 



Eq. 



for the action differences occuring in tnmt*, 



(with turn from Eq. (^)) and the sum rule (|^) one finds 
for this 'loop' contribution 
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with 



I{e,T)=ReJ\eP{e,T)e^p(^^^Y (11) 

In the semiclassical limit (ft — > 0) the contribution from 
small angles is dominant. In Eq. (pi]), the density P(e, T) 
of self-crossings with angle e for a long orbit of time T 
can be expressed as an integral over all loops, associated 
with the self-crossings, with times Tj^ini^) < t < T: 



dt{T - t) sin(e) p, 



erg , 



(12) 



where pcrg = l/(27rmA) is the ergodic classical return 
probability. The lower cutoff accounts for the fact that, 
owing to the exponential divergence of neighboring or- 
bits in a hyperbolic system, a minimum time Tmin(£) is 
required to form a closed loop from two trajectories start- 
ing at the crossing with initial angular difference e. Hence 
'7min(£) can be estimated from c ~ eexp(ATmin(£)/2) 
with c of order tt. Detailed numerical and analytical stud- 
ies |l^ have shown that this indeed holds true and that 
the number of crossings for T ^ cxd is given by 
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with Tmin(e) = -(2/A)log(e/c). The integral (|11]) over 
the leading-order term in Eq. (^ is purely imagi- 
nary, and thus its contribution vanishes. However, the 
contribution to /(e,T) of the second, logarithmic term 
in Eq. ( p^ ) is finite and gives —(h/2mA)T, independent 
of A. We then obtain from Eq. ( [lO|) 
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Hence, the lack of short loops with t < T,„in(£) gives rise 
to a negative quantum correction to the transmission. 

Correspondingly, we find for the loop correction to the 
reflection coefficient 



''nm lloop 



1 



{N + N'f 



(15) 



Here, as for the diagonal contribution (||), backscattering 
into the same channel is twice as probable. 

Summing over all initial and final channels we obtain 
for the leading-order quantum transmission and reflec- 
tion (5rioop = -NN'/{N + N'f and SRxoop = -N{N + 
1)I{N + N'f. For iV, iV' > 1 we have (5i?diag + -^-Rioop ^ 
NN' /{N + N'Y = -5T\oop- This imphes conservation of 
the average current in the semiclassical limit. Consider- 
ing off-diagonal terms allows us to semiclassically com- 
pute WL corrections consistently either in transmission 
or reflection. They precisely coincide with the RMT re- 
sult JTrmt = — 1/4 for iV = A^' ^ oo. Comparison with 
the RMT resuhs for finite N, N' suggests that the fcth 
order terms in Eqs. (|[||) correspond to semiclassical It- 
loop contributions; the diagonal terms are considered as 
0-loop and the orbits in Fig. |l| as 1-loop terms. 

Since the closed loops formed by the off-diagonal or- 
bit pairs are traversed in opposite directions, see Fig. [|, 
these orbits acquire an additional action or phase differ- 
ence in the presence of a weak magnetic field B due to 
the flux enclosed. For a uniform perpendicular field the 
action difference is AnAB/ipo, where A is the area of the 
loops and (po the flux quantum. We assume that the dis- 
tribution p{t; A) of enclosed areas for trajectories of time 
t is Gaussian with a system specific parameter /?, 
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This is usually well fulfilled for chaotic systems 
For finite _B-fields we have to perform an additional 
integration of the field-induced phase differences over 
the area distribution: J^^dAp{t; A) cos{4:Tt AB /(j)o) = 
exp(— t/is), with the magnetic time is = (47^/3i3/0o)^^• 
Up to timescales Tniin(£) a negligible flux is enclosed by 
loops with small crossing angles. We consider this by a 
respective time shift when inserting exp^—t/ts) into the 
integral (|l^) over loop lengths: 
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In Eq. ([l7| ) we used Tiiiin(£) ^ ^s- This corresponds to 
the original assumption, T'inin(£) ^ in the range of 
interest, t ^ ts- Only the term linear in rinin(£) con- 
tributes to the integral (|ll]), and we eventually obtain, 
after computing the T-integral (|l^), a Lorentzian field de- 
pendence of the transmission coefficient: |^nm(P)|^oop — 
|tnm(0)|ioop/(l -I- t/^b), a corresponding resuh applies 
to |''nm(-B)lk)op- This coincides with the Lorentzian WL 
lineshape obtained in the diagonal approximation 
making clear why the diagonal terms already qualita- 
tively account for the WL peak profile. The entire WL 
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correction from the diagonal and ofFdiagonal (1-loop) 
contribution then reads, in terms of the classical reflec- 
tion and transmission coefficients rd and ici, 
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Our refined semiclassical approach to the Landauer con- 
ductance yields the correct WL magnitude andlineshape. 

The Kubo conductivity reads, in terms of advanced 
(retarded) Green functions G^{r,r';E) (for a system 



of area A), a 



^n/(47rA)]Tr{z):rAGi4AG}, with 



AG = G+ — G^ . The trace is semiclassically evaluated in 
position representation by approximating the products of 
Green functions involved through double sums over clas- 
sical phase-carrying paths. Pairing identical orbits in the 
diagonal approximation leads to the classical Kubo con- 
ductivity [^j7|J|]; off-diagonal terms are again required 
to compute WL for chaotic systems. As a prototype of 
an extended clean chaotic system consider, e.g., a two- 
dimensional Lorentz gas. This has been experimentally 
realized by regular arrays or disordered ensembles of anti- 
dots in two-dimensional semiconductor heterostructures 
pl[ . The antidots act as classical scatterers giving rise to 
diffusive motion on long time scales, while the dynamics 
for intermediate times is governed by chaotic scattering. 

Our semiclassical treatment of WL is based on off- 
diagonal pairs of paths which have much in common with 
the orbits discussed above (Fig. 1(b)): they consist of one 
long self-intersecting trajectory being backscattered after 
multiple bounces with antidots with nearly opposite mo- 
mentum and a neighboring orbit which follows the loop 
formed by the first in opposite direction. A careful treat- 
ment of the conductivity trace integral for such pairs of 
backscattered paths (involving again cutoff times loga- 
rithmic in the crossing angle) gives rise to a non- vanishing 
negative quantum contribution 6(t at B — [p^ . 

This WL correction for chaotic systems with classical 
scatterers turns out to coincide with that from disordered 
systems with quantum impurity scattering. We find 
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(19) 



where t^ is the phase-coherence time and tc\ the elastic 
scattering time due to refiections at the antidots. Diffu- 
sive motion on long time scales, accounted for in a sum 
rule similar to Eq. (||), is reflected in the In. Eq. ( |l9| ) 
coincides with the result of Ref. [|j for antidot systems 
when tfj, is large compared to the Ehrenfest time. 

The approach above can be generalized to treat linear- 
response functions of other observables ||l^ . 

To conclude, a semiclassical treatment beyond the di- 
agonal approximation is appropriate to compute quan- 
tum corrections to the average conductance in clean 
chaotic conductors, both in the Landauer and Kubo 
framework. Chaotic classical dynamics is responsible 
for a logarithmic angular dependence of the classical re- 
turn probability, respectively the loops involved, which 



turns out to be crucial for computing weak localization. 
Numerical results for billiards |lj] show that this loge- 
dependence holds true also for nonuniformly hyperbolic 
systems indicating that the mechanism presented here is 
rather general. A semiclassical evaluation of higher-order 
loop corrections is yet to be performed. While such terms 
are not negligible for the spectral form factor, the one- 
loop corrections considered here play the dominant role 
for quantum transport in the mesoscopic regime. 
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